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Angular-resolved photoemission data on half-filled doped cuprate materials 
are compared with an exact- diagonalization analysis of the three-band Hub- 
bard model, which is extended to the infinite lattice by means of a perturbation 
in the intercluster hopping (cluster perturbation theory). A study of the band 
dispersion and spectral weight of the insulating cuprate Sr^CuO^C^ allows 
us to fix a consistent parameter set, which turns out to be appropriate at 
finite dopings as well. In the overdoped regime, our results for the spectral 
weight and for the Fermi surface give a good description of the experimental 
data on Bi2Sr2CaCu20$+s ■ In particular, the Fermi surface is hole-like and 
centered around k = (vr,7r). Finally, we introduce a hopping between two 
layers and address the issue of bilayer splitting. 

PACS numbers: 74.72.-h, 74.72.Hs, 79. 60. -i. 

1. INTRODUCTION 

Angle-resolved photoemission spectroscopy (ARPES) has provided sig- 
nificant insight into the single-particle properties of high-temperature su- 
perconductors (HTSC), providing a measure of the /c-dependent spectral 
function below the Fermi surface, which is suitable for a direct comparison 
with theoretical models. Recently, particular attention has been devoted to 
the spectral properties of undoped antiferromagnetic parent compounds of 
the HTSC, mainly Sr2Cu02Cl2 due to its ideal experimental characteristics. 
A series of photoemission measurements E&B of these half-filled cuprates pro- 
vides a detailed view of the quasiparticle dispersion of the highest electron- 
removal state, commonly known as the Zhang-Rice-Singlet (ZRS) @. In a 
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number of previous works, the wave- vector dependence of this band has been 
theoretically analyzed on basis of exact-diagonalization (ED) and quant 
monte-carlo (QMC) studies of Hubbard t-J B and three-band models Q' 
On the other hand, it is clear that an appropriate inclusion of the most 
relevant orbitals (Cu3d a; 2_j /2 , 02p x and 02p y ) of the Cu02 unit cell, i. e., 
based on the three-band Hubbard model should provide a more complete 
description of the system. The trouble is that an exact-diagonalization anal- 
ysis of this model has to be necessarily restricted to a rather limited number 
of unit cells, due to the rapid increase of the size of the Hilbert space. Even 
QMC simulations of larger systems were still severly restricted, providing 
informations on only a few k-points in the Brillouin zraie. 

To remedy this problem, we follow here the idea Ii3 of extending the re- 
sults of the exact diagonalization of a small cluster to the whole lattice within 
a lowest-order perturbation in the intercluster hopping. This technique has 
been shown to produce quite good results for the Hubbard 1121 model and 
has been successfully applied to the study of the stripe phase in the cuprate 
materials IHrta. It has the advantage to treat accurately short-range correla- 
tions, which are believed to be the most important ones in these systems 0. 
Furthermore, at the same time, it accesses all k points of the Brillouin zone, 
while in ordinary ED only few k points are available. This last property is, 
of course, crucial if one wants to analyze a dispersion or determine a Fermi 
surface with sufficient details. 

As a matter of fact, recent advances in experimental ARPES tech- 
niques allowed for a direct mapping of the Fermi surface, which has been 
analyzed in detail in overdoped Bi2Sr2CaCu20s+5 (Bi2212) 
l. Unfortunately, there seems to be disagreement between 
results obtained at different photon energies. Measurements taken at a pho- 
ton energy of 22eV typically show a Fermi surface closed around k = (ir, it) 
("hole-like"), whereas at 33eV photon energy the same materials show a 
shape closed around k = (0,0) ("electron- like"). To understand this is- 
sue, we extend our analysis to an overdoped system with 25 % doping, 
where we keep the same parameters obtained at half filling. Our results 
give a quasiparticle dispersion which is compatible with the experiments on 
the overdoped cuprates, exhibiting a hole-like Fermi surface closed around 
k = (it, 71"). 

In the final part of the paper, we present results obtained by coupling 
two layers with an inter layer hopping t± = 0.1t p ^, and study its effects on 
the Fermi-surface splitting, comparing with experimental results. 
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2. MODEL AND TECHNIQUE 

Common property of all HTSC and of their parent compounds is the 
presence of stacked Cu02 layers. The general agreement is that the rele- 
vant excitations of these materials take place within these layers. A generic 
model considers three relevant orbitals, namely, Cuj „ „ , 0„ and 0„ 
orbitalsO'EjlH3. Restricting to the on-site Coulomb interactions withinthe d 
copper orbitals only, this gives the three-band Hubbard Hamiltonian c3 

(y},<7 

-tpp a 'jj'{p],aPj'^+Pf,aPj,<r) 

+yE d lA,° d la> d i,«' 

+ A E^W' (!) 

where the operators d\ a and pj a create holes in the Cu 3d and in the O 2p 
orbitals, respectively, and and a'--, give the usual orbital phase factors. 
(...) denotes summation over nearest neighbours. 

Exact diagonalizations of are usually limited to 4 unit cells, allow- 
ing for just 3 inequivalent wave vectors, which are not enough to address 
questions of quasiparticle dispersion or of Fermi surface. Nevertheless, since 
the important correlations of these strongly-correlated materials should be 
relatively short range (a celebrated example is the Zhang- Rice singlet 0) one 
can still envisage to capture these short-range correlations by exact diago- 
nalization of a small cluster, and then continue the cluster properties to the 
infinite lattice within a perturbation in the inter-cluster hopping elements. 
This idea has been suggested by Senechal and coworkera!2l and has been 
termed "Cluster- Perturbation Theory" (CPT). A closely related method has 
been previously suggested to treat the coupling between one-dimensional 
Luttinger liquids, ancLtQ address the problem of the crossover from one to 
higher dimensions EScJ'Bj. 

Central piece of this method is the Green function of the cluster Gij a (z), 
which is obtained from Lanczos diagonalization. At the lowest order in the 
CPT, the full cluster Green's function is non-diagonal in the momenta, as 
the partition into clusters produces a smaller "superlattice" Brillouin zone 
(BZ). More specifically, one has 

1 cluster 

GS+G^+G 2 ,A Z ) = If £ e^-^e* G ^-* G ^g ri , r2 A<l,z), (2) 

ri,r 2 
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where Gj are vectors of the reciprocal superlattice, rj are cluster sites, q is 
restricted within the superlattice BZ, and Nq is the number of unit cells of the 
cluster. Here, the cluster Green's Function in the superlattice representation 
^ri,r 2 ,o-(q> z) can be written as 

£ ri , r2lCT (q, z) = iG^z)- 1 - V(q)j ~* , (3) 

L J ri,r 2 

where both the cluster Green's function G a (z) as well as the Fourier-transformed 
intercluster hopping V(q) are Nq x Nq matrices in the cluster-site indices. 
y(q) is the only term connecting different clusters and, thus, displaying a 
q-dependence. More specifically, let T riiI . 2 (A) be the amplitude of hopping 
of a particle from a site ri in cluster R and a site T2 in cluster R+ A (which, 
for a translation-invariant system, only depends on the relative position, A, 
of the two clusters) then V(q) is given by 

nq) ri , r2 =E e " q ' Ar r 1>r2 (A). 
A 

This method obviously allows for a perturbative treatment of intracluster 
hopping terms as well. These can be included via the intercluster part 
T rijI . 2 (A = 0). This gives the possibility of diagonalizing the cluster with 
periodic boundary conditions (BC), which is, clearly, more convenient nu- 
merically, but which contains unphysical hopping terms within the cluster. 
These additional terms, which are not present in the original problem, are 
then removed perturbatively via a T rijI . 2 (A = 0) term with the opposite 
sign. 

The non-diagonal terms Gi / G2 in Eq. (|2|) are produced by the dif- 
ferent treatment of the intra- and intercluster hopping terms. Within the 
CPT, these terms are neglected, and one takes only the Gi = G2 part of 
the Green's function to evaluate the spectral properties. Of course, in our 
multiband Hubbard model, the Green's function is additionally a 3 x 3 ma- 
trix in the orbitals of a single cell. The spectral function plotted in figures [l] 
and I is given by the sum of the diagonal elements of the Green's-function 
matrix. 

Notice, that the lowest-order CPT becomes exact in the trivial case of 
vanishing intercluster hopping and for non-interacting systems. This is due 
to the fact that corrections to Eq. (§) with Eq. @ are given by higher-order 
cumulants, which are vanishing when Wick's theorem holds, i. e. for non- 
interacting electrons. This fact makes the CPT an appealing interpolation 
between the strong- and the weak-coupling limits. 

Our cluster consists of Nq = 2x2 unit cells of the 3-band model, for 
a total of 12 lattice sites. For a single doped hole, the calculation of the 



Spectral Properties of High-T c Cuprates 

Green's Function requires 156 diagonalizations of with at most 48400 
basis states. 

In principle, if we require that the hole density of the cluster fixes the 
doping of the system, the lowest doping we that can treat (beyond half 
filling) is 5 = j. Even so, we have to average the Green's function between 
the two situations where the additional particle in the cluster has spin up or 
down. 

As a matter of fact, one can also treat any doping between 5 = and 
5 = 25%, by considering the cluster as being in a grand-canonical ensemble 
with fixed chemical potential [i between hole numbers = 4 and iih = 4+1. 
For finite temperature T, one can always fine-tune /i so that the average 
is the desired one = 4(1 + 5). In the T — > limit, this is obtained by using 
as a cluster Green's function to insert in Eq. (^) the weighted sum of the 
two Green's functions at = 4 and = 5 (G n =4 and G n= 5, respectively), 
i. e., G <nh>=i+X = x G nh=5 + (1 - x) G nh= 4 ( with < x < 1), at 
the appropriate chemical potential. This is, physically, what one expects 
to occur in a cluster embedded in a larger system with S = § , namely, the 
particle number of the cluster n will mainly fluctuate between n = 4 and 
n = 5. 



3. RESULTS 

First, we will discuss the dispersion of the highest electron removal 
states on the basis of the most recent experimental studies of the half-filled 
cuprate Sr2Cu02Cl2. Up to few exceptions i, these studies provide a consis- 
■Jit picture of the quasiparticle dispersion of the HTSC parent compounds 
3, showing two parabolas centered at k = (ir/2,0) and k = (tt/2,tt/2). A 
characteristic energy difference A a ~ 100 — 150meV between the maxima 
of the parabolas is observed. The total width W of the first band is about 
300 — AOOmeV (compare fiffure_0 for experimental data). In accordance with 
earlier numerical studies I3BBM1I of the three-band model, we find that the 
parameter set t pp = 0At p d , U = 6t p d and A = 4t p d gives the best agreement 
with the experimental band structure. With the value of the copper-copper 
hopping matrix elementEil taken from the literature, t p( i = 1.3eV, our cal- 
culations yield A a = 150meV and W = -iOOmeV. For the charge transfer 
gap we get At = 1.8eV in good agreement with the experimental value 
of 1.5 — 2eV Ej@. All these quantities are well within the margins of the 
experimental constraints. Therefore, these parameters should give a rather 
good description of the materials. Figure [l] shows a density plot of the 
single-particle spectral function A(k,uj) along the high-symmetry directions 
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Fig. 1. Grayscale plot of the spectral density of the half filled three-band 
Hubbard model, Eq. ([[]), with t pp = 0.4, Ud = 6 and A = 4 (in units 
of t p d = 1.3eV), as obtained from CPT (left). The points mark the peak 
positions. The right rjanel shows theband dispersion obtained from different 
experiments (circles a, filled circles and triangles □). 



of the Brillouin zone (left panel). On the right panel we show corresponding 
ARPES data for comparison!0B. The dispersion of the band is reproduced 
quite well by our calculation. On the other hand, the spectral weight is too 
strong in the region (tt,0) to (vr,7r), and too weak in the region between 
(0,0) to (it, 0). This apparent disagreement is discussed in more detail be- 
low. Note that the scale of the binding energy is not absolute, since we are 
dealing with insulators. 

We now comment on previous analyses of the experimental band disper- 
sion displayed in Fig. |l| which have been based on the t-J model. This model 
is believed to be an effective low-energy version of the three-band Hubbard 
modelil. The t-J model is more appropriate for numerical analysis, since his 
Hilbert space is smaller for a given number of unit cells. However, our results 
obtained by a CPT analysis of the three-band Hubbard model clearly differs 
from what one generally obtains from the t-J model. The latter predicts, 
in contrast to experiments, a flat band close to the Fermi energy at (ir, 0) 
and (tt, 0) — (0, 7r) 0. A solution to this shortcoming of the t-J model has 
been suggested, by introducing longer-range hopping and interaction terms. 
In fact, a mapping of the three-ban dHubbard model into a generalized t- 
J model has been recently proposed EJ, which is consistent with results of 




f-consistent Born approximation calculations and experimental data 
. In particular, these latter results exhibit a dispersion in excellent 
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Fig. 2. Spectral density of the overdoped (5 = 0.25) three-band Hubbard 
model as obtained by CPT (left panel). Parameters are as in Fig. |. The 
white dots mark the peak positions. Corresponding experimental results 
from Ref . ^ are included for comparison (right panel) . Empty circles denote 
the overdoped, filled circles the underdoped case. The solid line marks the 
CPT result. 



agreement with our calculations of Fig. |]. 

From Fig. |], one can see that our calculation underestimates the spectral 
weight along (0,0) — (vr,0), which is usually strong in ARPES. One reason 
could be related to the small cluster in which the self-energy is evaluated, 
which takes only partially into account the electron correlations. Moreover, 
it was established numerically 113 that quasiparticles contained in the cluster 
are reproduced correctly, but they loose coherence over the cell boundary. 
Alternatively, the different spectral weight could be caused by intrinsic limi- 
tation of the model. For example, neglected orbitals " - although distant in 
energy - could possible give significant contributions to the spectral weight. 

We now turn to the doped systems. Here, we assume, as a first guess, 
that the parameters are not significantly dependent on the doping and on the 
particular material, and we keep the same ones for the doped region. The 
minimum doping that can be achieved with a 12-site cluster by adding an 
additional electron is 5 = 0.25. This undoubtedly high doping is nevertheless 
well suited to be compared to the detailed experimental results on dispersion 
E and Fermi surface 000000 of the overdoped cuprates. The results 
of our calculation are plotted in Fig. ||. The quasiparticle dispersion from 
k = (0, 0) to k = (-7T, n) exhibits a steep ascent beginning at a binding energy 
of ~ — 0.5eV and crossing the Fermi energy shortly before k = (tt/2,tt/2). 
Starting from k = (0,0) and going along (1,0), the dispersion reaches a 
saddle point at ~ k = (vr,0) well below the Fermi energy at ~ —WOmeV. 
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Fig. 3. Spectral weight at the Fermi energy giving an estimate of the position 
of the Fermi surface. Results are displayed for a single Cu02 layer (a) and 
for the bonding (b) and antibonding (c) band of a bilayer with interlayer 
hopping t± = O.ltpd- The white dots indicate the maxima of the spectral 
weight on a given branch. Panel (d) shows the comparison with experimental 
results on overdoped (Bi,Pb)-2212 at 22eV photon energy H(empty squares) 
and and slightly overdoped (Bi,Pb)-2212 at 55eV photon energy H3(empty 
circles) . Filled squares and circles denote the results of the CPT calculation 
for the bonding and antibonding branch, respectively. 
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A similar saddle-point dispersion is observed in the experimental data for 
Bi2Sr2CaCu20g+5 0, shown for comparison in the right panel of Fig. The 
good qualitative agreement of our numerical results with the experimental 
ones suggests that our approach of taking the same parameter values in the 
undoped and in the overdoped systems is Quite reasonable. This was also 
found earlier by intensive QMC simulation a. 

The magnitude of the spectral weight at the Fermi energy is displayed 
as a grayscale plot in Fig. ||| for the whole Brillouin zone. The darkest 
regions, thus, give an estimate of the location of the Fermi surface. The 
data clearly show a hole- like barrel closed around k = (tt. n)^as usually 
observed in ARPES experiments on the overdoped cuprates I15E3l3E3, mostly 
at a photon energy of Eu,~ 22eV . On the other hand, experiments carried 
out at a higher energy HtlEil E v = ?>2>eV disagree with these results, and 
report an electron-like Fermi surface. This is quite puzzling, and there is not 
yet a full consensus about this issue. Some recent papers suggest that the 
observation of two different Fermi surfaces is actually due to the splitting 
from the coupling between the bilayers EjU. Nevertheless, our data suggest 
that, if the bilayer splitting is not relevant, the Fermi Surface should be 
hole-like. 

The coupling between two layers can be treated without further diffi- 
culties within our method. Again, we diagonalize a 2 x 2 unit cell exactly 
and include an interlayer hopping term t± within the CPT treatment. Since 
the effective hopping t± under consideration is relatively small and CPT 
amounts to a perturbation in the intercluster hopping, we believe that this 
is a good procedure. 

The interplane hopping has been analized in studies based on LDA cal- 
culations and has been mapped on effective multi-orbital models by integrat- 
ing out high-energy degrees of freedom E3. This approach provides results 
which are consistent with the experimentally observed bilayer splitting El. 
In our calculation, we adopt this idea by treating the orbital phase factors of 
the various interplane hoppings accordingly. We take a value of t± pa 0.1t p d, 
which is sufficient to produce an observable splitting of the Fermi surface. 

Figure || (right panel) displays our results obtained for the bilayer. Dif- 
ferent experimental results are shown in the last panel for comparison. Two 
branches can be distinguished, one being almost a square closed around 
k = (w, 7i"), The other one is closing around k = (0,0), forming an electron- 



like Fermi surface as observed in Ref. 19. Similar results have also been 
found by first-principle calculations of ARPES in Bi2212 El. These results 
suggests, that a possible reason for the crucial differences in the experimen- 
tal observations might be found in interlayer processes. An explanation of 
the apparent dependence of the Fermi surface on the photon energy would 
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require the consideration of matrix element effects which have not been ad- 
dressed within this paper, but are treated in detail elsewhere H. 

4. CONCLUSION 

This paper reports a systematic study of single-particle spectral prop- 
erties of the three-band Hubbard model in the half-filled and in the doped 
region. In order to obtain sufficiently accurate results and to sort out details 
of the band dispersion we adopt a method (so-called CPT) combining ex- 
act diagonalization of small clusters with a perturbation in the intercluster 
hopping. This analysis allows us to determine a parameter set, which appro- 
priately describes both the insulating as well as the overdoped region and 
gives a good agreement with ARPES data. The Fermi surface obtained at 
high dopings is hole like, although an interlayer hopping of about t± sa 0.1t p d 
produces a splitting into an electron and a hole-like Fermi surface , in agree- 
ment with ARPES experiments. 

ACKNOWLEDGMENTS 

We thank R. Eder for useful discussions. This work was supported 
by the projects DFG HA1537/17-1 and HA1537/20-1, KONWIHR OOPCV 
and BMBF 05SB8WWA1. Support by computational facilities of HLRS 
Stuttgart and LRZ Munich is acknowledged. 

REFERENCES 

1. C. Diirr et al, Phys. Rev. B 63, 014505 (2001). 

2. S. LaRosa et al, Phys. Rev. B 56, R525 (1997). 

3. B. O. Wells et al, Phys. Rev. Lett. 74, 964 (1995). 

4. F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 (1988). 

5. R. Preuss, W. Hanke, and W. von der Linden, Phys. Rev. Lett. 75, 1344 (1995). 

6. R. Eder, Y. Ohta, and G. A. Sawatzky, Phys. Rev. B 55, R3414 (1997). 

7. G. Dopf, A. Muramatsu, and W. Hanke, Phys. Rev. B 41, 9264 (1990). 

8. G. Dopf et al, Phys. Rev. Lett. 68, 353 (1992). 

9. G. Dopf et al, Phys. Rev. Lett. 68, 2082 (1992). 

10. D. Senechal, D. Perez, and M. Pioro-Ladriere, Phys. Rev. Lett. 84, 522 (2000). 

11. M. G. Zacher et al, Phys. Rev. Lett. 85, 2585 (2000). 

12. M. G. Zacher et al, cond-mat/0103030. 

13. P. Fulde, Electron Correlations in Molecules and Solids (Springer, Berlin, 1995). 

14. H. Fretwell et al, Phys. Rev. Lett. 84, 4449 (2000). 

15. A. D. Gromko et al, cond-mat/0003017 (2000). 



Spectral Properties of High-T c Cuprates 



16. D. L. Feng et al, Phys. Rev. Lett. 86, 5550 (2001). 

17. J. Mesot et al, Phys. Rev. B 63, 224516 (2001). 

18. M. S. Golden et al, Physica C 341 348, 2099 (2000). 

19. P. V. Bogdanov et al, cond-mat/0005394 (2000). 

20. D. S. Marshall et al, Phys. Rev. Lett. 76, 4841 (1996). 

21. W. Brcnig, Phys. Rep. 251, 153 (1995). 

22. Z.-X. Shcn, and D. S. Dessau, Phys. Rep. 253, 1 (1995). 

23. V. J. Emery, Phys. Rev. Lett. 58, 2794 (1987); V. J. Emery and G. Reitcr Phys. 
Rev. B 38, 4547 (1988). 

24. E. Arrigoni, Phys. Rev. Lett. 80, 790 (1998). 

25. D. Boies, C. Bourbonnais, and A.-M. S. Tremblay, Phys. Rev. Lett. 74, 968 
(1995). 

26. E. Arrigoni, Phys. Rev. Lett. 83, 128 (1999). 

27. S. Uchida et al, Phys. Rev. B 43, 7942 (1991). 

28. S. L. Cooper et al, Phys. Rev. B 47, 8233 (1993). 

29. M. S. Hybertsen et al, Phys. Rev. B 41, 11068 (1990). 

30. T. Tohyama and S. Maekawa, Supercond. Sci. Technol. 13, R17 (2000). 

31. J. M. Eroles, C. D. Batista, and A. A. Aligia, cond-mat/9812325 (1998). 

32. F. Lema and A. A. Aligia Physica C 307, 307 (1998). 

33. A. Damascelli, D. H. Lu, and Z. X. Shcn, J. Electron Spectr. Relat. Phenom. 
117-118, 165 (2001). 

34. O. K. Andersen et al, J. Phys. Chem. Solids 56, 1573 (1995). 

35. A. A. Kordyuk et al, cond-mat/0104294 (2001). 

36. S. Legner et al, Phys. Rev. B 62,154 (2000). 

37. Y. D. Chuang et al, Phys. Rev. Lett. 83, 3717 (1999). 

38. A. Bansil and M. Lindroos, Phys. Rev. Lett. 83, 5154 (1999). 

39. C. Dahnken and R. Edcr, cond-mat/0109036. 



